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SYNOPSIS 
(Part a) 

In this part of the worki the analog simulation 
of a hydro-electric governor system is presented, taking 
i.'ito account the true turbine characteristics. Effect 
of parameter variations such as speed dependence of flow, 
gate variations are studied both on stability and 
transient response of the system, and optimum adjustment 
of governor parameters is indicated. Comparison is 
made between the different models of the system including 


the classical model 
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SYI^OPSIS 
(Part B ) 

In this part of the work the m^dal control approach 
is empl»^yec? to improve the transient response of a power 
system hy shiftin^i the dominant eigenvalues of the un- 
controlled system farther towards the left in the complex 
plane, using complete as well as incomplete state feedback. 

A sector criterion is used for preassigning the new 
locations of the eio^nvalues and selection of appropriate 
control inputs for shifting is indicated. The advantages ■ ^ 
of the modal control approach over the frequently used 
approach using optimal linear stote regulator formulation ‘ 
are discussed.' 
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CHAP TEE 1 


BTTRODUCTIOH 


The study of optimum adjustment of governor paraLictcrs 
of a hydro-electric system, hitherto, have used in most 
cases an approximate model for the turbine part of the 
porer system. The true turbine characteristics were 
completely ignored, which as shown by Goldwag [ij have 
significant effect on the stability and transient response 
of the system. Step by step methods are available, which 
take into consideration the true turbine characteristics 
and give solutions with greater accuracy. But these 
methods are too involved when dealing with problems with 
small excursions of gate or problems of stability. Hence 
a labour saving method, representing at the same time the 
turbine more accurately, is needed to have better results. 

A more accurate linearized model of water turbine 
has been developed by Goldwag [l] and theori tical investi- 
gation (Appendix B) of the problem of stability for all 
types of machines at all gate openings has been made by 
him. He has also discussed the effect of turbine 
cnaracteristics on the optimum setting of governor 
parameters based on two criteria. According to him, if 
we draw the limit of stability curve, and if we denote 
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the point of critical damping and further draw the curve 
corresponding to the product (temporary droop and 

dash pot time constant product), obtainable at the point 
of critical damping, then the possible region of operation 
would be between the two curves. The former cirrve does 
.not represent a preferred setting of governor parameters 
and the latter critical damping curve may or may not 
represent a suitable setting dt;pending on the value of 
(hj-T, ) . In any case, a setting of governor parameters 

0 a Qr 

wj th b,T, value less than or equal to (b . T, ) is 

t d T d. QP 

nox-ma.lly preferred. 

In a recent paper Pai and Gfanesan [2] have 
discussed the optimum adjustment of governor setting 
using parameter plane technique. They have used only an 
apx'roximate model for turbine. It will bo shown in the 
Chapter 2 that their model which is referred as the 
classical model is a special case of that of Goldwag. 
Purther more compatibility between these models is 
established via the well known block diagram language 
of the control engineer. 

Another accurate linearized model of the turbine 
is given by Ramey and Skooglund [ 3 ]. This suffers 
from the dffect that certain coefficients are to be 
determined experimentally for the particular 
installation under consideration. In Chapter 2 a 
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comparison is made between this model and that given by 
Goldwag. 

In Chapter 3 analog simulation of the hydro-electric 
governor system, as given by Goldwag is carried out and the 
theoretical results of stability as predicted by Goldwag 
an: verified. In addition the effect cf self-regulation 
of the system, which was neglected by Goldwag in his 
theoretical investigation, has been studied and the results 
•verify with those given by Pai and Ganesan [2]. 



CHilPTLE 2 


HIDHO^liLLCTRIC G-OlhimOR SYSTEM^ MODPL 

2,1 Math. epiati cal Model takln,^ into account^ True lur'bin.e 
Characteristics 

The basic differential equations of the water 
turbine system as such can be derived from the information 
of the turbine characteristics like rela,tive unit quantity 
(RUT) versus relative unit speed (RbS), relative unit 
torque (RUT) v- rsus RUu and relative servo opening (RSO) 
versus RUQ and ROT. These characteristics are available 
either as experimental dS'ta or can bo derived from Mussel 
Curves [l,4]» usually provided. The basis for the 
derivation of the turbine S 3 ?"stem differential equations is 
that the above mentioned characteristics can very closely 
be represented by straightline tangents at the point 
corresponding to the unit speed around which an excursion 
of gate or the stability is to be investigated . As a result 
certain new parameters are introduced (Figures A.3» A.4» 
A.5)j which are 

i) P ~ the slope of RUQ versus RUS characteristic 

ii) ~ slope of RUT versus RUS characteristic 

iii) a - the relative effective gate, which is defined 


by the equation 

a = g + Y s 


( 1 ) 
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vhero y is the slope of effective gate versus RSO 
cberactcristic (Figure A.5)* 

The detailed doriva.tion of the turbine system 
dynamic equations is given in Appendix A. The two 
hai-sic differential equations derived for the turbine 
system ore as follows 


w 


ind 


T 


dn 
m dt 


do 

dt 


= Y A s + 


j. (o^-- g) 


dt 


dt 


+ P 


dn 

dt 


Z - [(a - P) - 


eg]n 


(2) 

(5) 


svbject to the assumptions that the slope of the RUQ 
versus RUS characteristic does not vary significantly 
with gate and the bulk water hammer represents adequately 
the pressure changes following the variations of flow. 

Since the main interest is to show the effect of 
the true turbine chara,ctoristics , the most basic form of 
transfer f'unction for the governor is assumed and is 
giver by 


As 


bfTdS 


(4) 


Equations (2), (5) and (4) represent mathematically the 
isolated hydro-electric system under consideration. 
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2 . 2 Bloclr Diagram Represeiatation of th.e_ System 

The mathematical model of the system given hy 
oqns.(2), (5) and (4) can conveniently be represented by 
means of a block diagram, well known to control engineers, 
the development of which follows. 

The governor transfer function given by eqn. (4) 
can be represented by a block diagram given in Figure la. 
Tne corresponding state variable diagram is given in 
Figure lb. 

The turbine system eqii8.(2) and (3) can be 
rGX)resented by the two block diagrams given in Figures 
2a and 5a. The corresponding state variable diagrams 
are given in Figures 2b and J'b. 

Combining the blocks given in Figures la, 2a and 
5a and doing certain simple maniptilations , a composite 
block diagram given in Figure 4a is obtained, which 
represents the hydro-electric governor system. The 
corresponding state variable diagram is given in 
Figure 4b. In the Figures 4a and 4b, 6 is defined as 
effective temporary droop, which is given by 



and A 1 is a step load change. 



la; Hiock niagrrirn 


(b) State varinbU' diagram 


nr. 1 GOVERNOR RE PRE SENTATION 



(a) Block diagram (b) State variable diagram 

FIG.3 SYSTEM REPRESENTATION 
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(t;) St ah varialdc diagram 
F--IG-4 MYDRO-Fl.ECTrdC GOVERNOR SYSTEM 


Tm S+D 


CLASSICAL MODEL OF HYDRO-ELECTRIC 
GOVERNOR SYSTEM, , „v T 
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2 . 3 Oomparison_ vith Classical Model 

The conventional hlock diagram representation that 
considers approximate turbine transfer function (classical 
model of turbine), as used by Pai and Ganesan [2] is 
given in Figure 5 where D is considered as load damping 
factor. 

When the conventional block diagram is compared 
with that given in Figure 4a, which takes into account 
true turbine characteristics, the following conclusions 
confirm the identities. 

(i) At rated load the value of a = 1.0 (as per the 
Figure A.5). While approximate turbine model is considered, 
it is assumed that flow doos not depend on speed and so 

the value of p = 0, With these values of a and p, it 
can be seen that the true transfer function of turbine 
in Figure 4a will reduce to the classical turbine transfer 
function shown in Figure 5- 

(ii) Referring to Figure A,5» if if wore to be assumed 
as an approximation for turbine characteristics that the 
rate of change of RUQ and RUT with scrvo->positions is 
uniform throughout the range of servo- openings (instead 

of servo-positions around the region of peak efficiency), 
then the slope y of effective gate will be equal to one. 

In that case there is no need for defining effective 
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temporary droop 6 and tiie transfer function of governor in 
Figure 4a will be same as that given in Figure 5» 

(iii) In Figure 5» D is considered as the load damping 
factor, where as in Figure 4a, corresponding to D in the 
transfer function of the system we have e^ ^ the overall 
coefficient of self-regulation of the system compromising 
s elf -regulation of the turbine as well as that of load 
and is given by 

n t g ^ ^ 

2.4 Comparison with Another Model [ 3 ] 

^lother linearized accurate model of a water 
turbine is given by Ramey and Skooglund [5j» Tko 
equations of quantity and torque that define the' turbino- 
penstock transfer function are 


<1 = Z 

+ a^2 ^ ^ ^ 

(7) 

A m = a 22 _ Z 

+ ag2 I' + ajj Ab 

(8) 


where partial derivative of turbine flow with 

rosnect to turbine head, 

a ^2 partial derivative of turbine flow with 

respect to turbine speed, 

a^^ is the partial derivative of turbine flow with 
respect to gate position, 

a^^ is the partial derivative of turbine torque 


with respect to turbine head, 
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^22 partial derivative of turbine torque 

with respect to turbine speed, 

is the partial derivative of ttirbine torque 
with respect to gate position* 

The corresponding equations of quantity and torque 
that arc given by Goldwag [l] are 

q = Y^s + z + pn (9) 

Am = Y A s + Z - (a - p)n (10) 

It can bo seen that the two sets of eqns.(7), (8) 
and (9), (10) arc identical for the following relationships 
of their coefficients* 


H 

1! 

^23 ^ 

^•11 

( a - P )/2 

a 2 i = 

(3a-p)/2 

^12 

P 

^22 = 

= -0 


( 11 ) 


However, the equations given by Goldwag are more 
realistic in the sense that the coefficients are 
c; efined in terms of the turbine characteristics that 
will normally be available, whereas the coefficients in 
the eqns. (7) and (8) are to be determined experimentally 
for the particular installation tmder study. 


CHAPTER 3 


,SIFiULATIOI'' OF TH:6 HrPRO-RLJiiCTRIC GOVERNOR SYSTEM 

3 *1 RanM'e of Variations of Bifft.rent Parameters 

In order to simulate the hydro-electric governor 
system as represented in Pigure 4? the range of variations 
of different parameters aro to bo known, which are as 
foil ows ; 

(i) p - the slope of RUQ versus RUS characteristic: 
For low head water turbine machines the value of p 

ic grea.ter th^n zero and for high head machines p is Icsss 
than zero. Hence to cover all typos of turbines the range 
of p thj'vt can be considered accordingly is 

- 1.0 ^ p < + 1.0 ( 12 ) 

at; rated load when a = 1.0. 

(ii) a - the relative effective ga.te: 

The value of effective gate at rated load is equal 
to one. Since our investigation is aromd the gate at 
which peak efficiency is obtained, effective gr-te values 
around one aro considered. The range of values considered 
is from 0.9 to 1,05 as obtained from Pigure A.5- 

(iii) T^, T^ - Time constants of water and turbine 
respectively: 

Following typical values are assumed for T^ and 
T^ as considered by Pai and Ganesan [2] 
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- water starting time = 1,24 secs. 

- Machine starting time = 9.05 secs. 

The values of T„ and T^ depend on the particular 
installation for which stability analysis has to be 
ixtad c * 

(iv) e^ - the coefficient of self-regulationt 

Tho effect of can usually be neglected since 
the solf-regulntion is more system than characteristic 
d ..pendent and depending on the type of load carried the 
overall effect of turbine and load s elf '-reg'ulat Ions may 
become zero. The value of can vary from zero to ono., 
Gold'wag has neglected in his theoretical inv evStigation 
of stability of the' system, 

(v) 6, T^ - effective temporary droop and dashpot 
time constants: 

Particular values of 6 and T^ constitute the 
optimum setting of governor based on the criterion used 
for the same. To obtain the limit of stability curves, 

6 and T^ are to be varied simultaneously over a range of 
values. Por this a range from 0.1 to 2.0 for 6 and a 
range from 0 to 1.0 for T^ can be considered, T^ value 
can go up to 2 or 5 in. case the damping factor or 
self-r€'gulation coefficient is not neglected. These 
ranges of 6 and T^ are assumed based on the theoretical 
results as obtained by Goldwag[l] and Pai and Ganesan[2]. 



3 1 2 Analog Simulation Results 


Analog simulation of tiie hydro-electric governor 
system as represented in Figure 4 is carried out using 
AC-20 analog compute- r. The experimental set up is shown 
ir Figure 6. A step load change ^ 1 of 10 percent is 
considered , 


Theoretically the limit of stability curves, as 
derived by Goldwag, are defined by a relation (Appendix B) 
given by 


X 


2 



whore 




m 


aT 


-IT P K 
m 


T 



and 

K = Up; , - P) 
2a 


(13) 

(14) 

(15) 


(16) 


Numerical data derived from eqin,(15) is given in Table 1 
from which theoretical limit of stability curves for 
various values of p can be dra>m as shown in Figure 7* 


Critical damping conditions are derived by Goldwag 
(Appendix B), which are as follows 

9{-K-l)X^=^ (1-X^) (I-X 2 ) . (17) 

(1 - X-, )* 

5(2-1) (I-Z 2 ) « 5 ^ 

a ^ p 


( 18 ) 




Tabic 1; Data for Theoretical Liisit of Stability Curves at Rated Load 
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liquation (17) 'will define the curves corresponding to the 
critical damping points and eqns.(l7) and (18) together 
vhen solved, will give the governor setting values for 
critical damping. The governor settings corresponding to 
critical damping condition for three types of turbines, 
i.e., for three values of p at rated load when a = 1.0 
are calculated which are presented in Table 2. numerical 
data as calculated from eqn.(17) is also presented in 
Table 3 from which the critical damping curves can be 
drawn as shown in Figure 7. According to G-oldwag the 
optimum Slitting of govi^rnor for a particular turbine will 
lie some whore in between the limit of stability curve and 
the curve correspond ing to the critical damping point. 


Table 2: Governor Settings for Critical Damping Condition 


Typo of 
turbine 

Effective 
temporary droop 

6 

Inverse 

time 

of dashpot 

constant 

l/Tfl 

0.0 

0.4050 


0.1035 

-.0.5 

0.6370 


0.0723 

+0.5 

0.2600 


0.1610 


During the course of simulation of the system 
on analog computer, the limit of stability curves for 
various machines are obtained and also the effect of not 
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neglecting as well as the effect of changing a which 
means operating the machine at a load other than rated 
full load, are studied. The numerical data for limit of 
stability curves for three values of p at rated full load 
neglecting as well as for o^ = 1.0, is given in Table 4. 
Variation of a is also considered the numerical data of 
which is given in Table 5, Various graphs are plotted 
as shown in Figures 7, 8 and 9- Transient response 
curves for the frequency state variable of the system 
for typical settings of governor are given in Figures 10 
to 15- 

3 • 3 Comments on the Analog Simulation Results 

(i) It is observed that there is significant effect 
of the true turbine characteristics (p) on the region of 
stable operation of the turbine as illustrated in Figure 7. 
The limit of stability curves at rated load obtained by 
analog simulation of the system for three values of p 
agree with those obtained theoretically ^7 Goldwag.. 

(ii) Tho optimum setting of the governor parameters 
6 and Tj tased on any criterion will normally Da In tho 
region between the limit of stability curve- and that 
corresponding to the point of critical damping. 

(iii) It is obvious frpm the transient response 
curves given in Figures 10 to 13, that a particular 
governor setting on different types of turbines will 
give rise to different types of transient responses. 
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Table 5: Ixperincntal Data for Limit of Stabilix;/ Curves 

Loads Different from rotrd losd. 
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Relative servomotor opening 

F IG- A-5 RELATIONSHIP BETWEEN GATEJORQUE, 
AND FLOW at CONSTANT UNIt'sPEED 



— Limit of stability curve 

Theoretical curve corres- 
ponding for critical damping 

o Critical damping point 
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(iv) Por gate oponings different from that at rated 
lonrl, wo observe from Pigure 9 that the region of stable 
op(,r‘ition becomes smaller for hight^r values of gate 
ov.cning.: and becomes larger for lower values. 

(v) If the damping coefficient e^ was not to be 
nogluctod, we observe from Figure 8 that the region of 
stable operction increases. This result agrees with 
tnat obtained by Pai and G-anosan [2]. 



CHAPM 4 


COHCLUSION 


Dlaf^rammatic representation of hyc^ro-eloctric 
,t,;overrjor syetora, taking true turbine characteristics 
into .'Account, is developed making use of the mathe- 
ticticnl model given by G-oldwag [l] . Analog simulation 
stuf'ics arc carried out for a particular case using 
AC*-20 analog computer. The theoretical results obtained 
earlier by Goldwag [l] and Pai and Ganesan [2] are 
verified, showing the significant cfiect of turbine 
characteristics on the stability of the system,. 



APPiiKDIX A 


P^RIVATIOI:^ OF TUBBIPA SYSTPH TRAKSFM PLHCTIOIV 
M AKINU USA OF TURBII TE CHARACTERISTICS 

The 130,810 o’ iff crontial equations of the water 
t'^rbine systcrj art; derived by Golrlvag [l] from the 
ij-’f ornrition of tbo following turbine characteristics. 

(i) Kelative unit quantity (fiUQ) versus rel'"tivc unit 
speed (RUS). 

(ii) Hc.'lative unit torque (RUT) versus RUS 
(lii)R(.lative' ec-.rvo opening (RSC) versus RUQ and RUT 

v.'hich arc given in Figures A.l, A. 2 and A^5. It can be 
obb'trve.l that these turbine ch'^racteristics can very 
clobely be reprt sented by straight lino tangents at the 
point corresponding to the unit speed as shown in 
Figures A*3» A. 4 and A. 5* The derivation of the turbin 
system equ-'^tions are as follows. 

Referring to Figaro A. 3 the RUQ (q^) at given 
gate and RUS(n^) is given by 

<11 -t P (n^-1) (A.l) ■ 

and the rolrtivc quantity (q) at head deviation Z is 


q = q^ fl+Z 

= qj,^ Vl+Z + P(l+n - iTl+Z) 


(A. 2) 
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sulDetituting for i^ivcn by 
= (l+n)/fl+Z (a. 3) 

Tills bu rcarrangoc! as 

q = (9lp-P) Vl+Z + (A. 4) 

The conventional form of the oqii,(A,4) is 

Comparing tqns. (A. 4) and (A*5)> wc sec that 
•.qn,(A.5) completely neglects the speed dependence of 
fi-ov,'. Secondly, when oqn.(xl.5) is used to represent 
propellers and low head machines, it over-estimates the 
effect of pressure changes of flow. . When used to repr>.sent 
high hufd machines, it undor-estirnates the dependence of 
T low on the pressure. iiqn.(A.5) was first derived for 
Pt Iton machines and. hence it represents correctly the 
flow for such medium head machines. . 

The value of p - tho slope of RUQ versus RUS 
oil. 'ract eristic, represents to some scale a measure of 
the error in the dependence of the queantity on pressure 
changes and also represents to scale the speed dependence 
of flow. 

ha uation of Torque ; 

Referring to Figure A. 4 the RUT (m^) at a given 
pete and RUS is given by 
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~ (A, 6) 

vhcru - tho slope, roprosonts the coefficient of self- 
x“'.i(j;Ulr tion of the turhinc and is given hy 

^t ~ (a. 7) 

''iiorc R is the relative speed dovi?^tion at tho extrapolated 
nniaway speed and rated load. 

The relative torqne at head deviation Z is given hy 

m = m^(l+Z) (A. 8) 

= ^}:+^ ] 1£ nfl+Z (A. 9) 

The first term in eqn.(A,9) represents the dependence 
of torque on pressure and the second term the inherent 
self-regulation of the turbine. 

When relative presvsure change is small and if its 
effect on turbine self-regulation is neglected, eqn.(A«9) 
can be reduced to a more manageable form 


^1 — ^Ir (^d'KZ) —• 


vherc 


K = 1 + 


1 - 

2R 


(A. 10) 


Por Pelt on and medium head machines the value of 
K = 5/2 and therefore, the conventional form of torque 
equation is 
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~ ^Ir (a, 11) 

and for other machines the value of K is to he properly 
chosen. 

Relative effective gate 

hov ve introduce the concept of relative effective 
^^ato Fhich facilitates in conoecting the parameters of the 
derived eqns.(l.4) and (A. 10) to the gate, thereby 
rodacin/'; the number of parameters. 

For any turbine, as shown in Figure A. 5, it is 
seen that over a range of servomotor positions around 
the region of peak efficiency, the RUQ and RUT are 
approximately equal and have equal rate of change with 
servo position. How, if we define relative effective 
gate in this region by the expression 

a = g + YS (A. 12) 


where y is the slope of effective gate versus ROO, then 
we have 


“ = ^Ir = “ir = S + Ys 


and 


Y = 


A a 
A s 


iilr ^ ^^Ir 

AS AS 


(A.13) 


(A. 14) 



Relative 







■fT+7 

K't'lutive unit speed 

FIO. A-3 DERIVATION OF EQUATION OF FLOW 



FIG. A-4 DERIVATION OF TORQUE EQUATION 
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Pi i‘f er&ntial oquetions of the tur'bine system 

If the changes in quantity and torque are small, 
then the efficiency of the unit can he assumed to he 
constant throughout the changes and the torque and quantity 
chn.re ct eristics can he related hy 

H = ^ n (A. 15) 


from vhich the diffsrence equation 
Am = y 4 s + 2 - (a-j3)n 


(A. 16) 


can ho derived. Comparing eqns.(A.lO) and (A. 16) and 
making use of eqns,(A.15) and (A. 14), we get 


R = a/(a~j3 ) 


= a-P 

K = 1 . ^1=21 


( ) 
2a 


(A. 17) 
(A. 18) 

(A. 19) 


V/ith the help of eqns. (A. .4), (A, 10), (A. 14) and 
(A.16), ^^e can write the two basic differential equations 
of the turhino system as follows; 


2 dq d( A s) (a-p ) ^2 dn 

1“ = dt = ^ “dt— + ^ p dt 

¥ 


(A. 20 ) 


T 


m 


dn 

dt 


= YA s + z - C(a-p)-'eg]n 


(A. 21) 


subject to the assumptions that the slope of RUQ versus 
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KUS characteristics does not vary significantly with 
g;ito end tho hulk water hammer represents adequately the 
pressure changes following variations of flow. 

Differontial equation of governor 

The transfer function of a sensitive governor 
(infinite gain is assumed) equipped with temporary and 
pt.rmanent droops is 

1+T. S 

As = ^ (A. 22) 

SJ 

P 

b^+h 

where (-r — ‘^) T. is termed as *on load' time constant 
P 

and S is the Laplace oeorator. 

Since our main interest is to show the effect of 
the trucj turbine characteristics, wo assinnc the most 
basic form of transfer function of the governor (permanent 
droop neglected) as 

(1 + T, S) 

As = n (A. 23) 

The eqns.(A.20), (A. 21) and (A. 23) represent 
mathematically the isolated hydro-electric governor system 
under consideration. 
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APPMIDIX B 

Tf-iJi.OKi.TICia RESULTS AS ODTAIPED BY GOLBWAQ fl] 

LiiTiit of stability 

Oomloijiing eqns . (x. 20) , (A. 21) and (A.23)> ‘th-© 
chnxr. cteristic equation of the system can he derived as 

cx T 

(l-6K;i3-6(K-l)e^)] + 

(1 + 60 ^ - ) S + 4- = 0 (B.l) 

^ ^d ^d 

Coefficient of self-regulation, can normally be 
neglcctfjd. Then the system characteristic equation is 

6V. " 

(1 - L = 0 (B.2) 

d d 

For a third order system, the characteristic 
equation of which can be represented as 

+ Aj + A^ S + i'Q = 0 (B.3) 

the system can be stable if the following conditions 
(Eouth-Hurwitz) are satisfied. 

A^ > 0, 0 , A^2i 0, Aq> 0 and 




(B.4) 
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iipj)lyif'i.^ thoae conditions to the system characteristic 
c lUfition (E,2) and observing tha.t at the limit of stability 
boeindory ^'-o' expression, which defines the 

limit oi stability curve for the system, can be derived as 


X 


2 


1 - X^ 


(B.5) 


1,'here 

a T a T 
> m m 

a T 

_ (B.7) 


end K is as given by eqn.(A.19). The expression (B.5) 
di-fines the limit of stability curves for various values 
of p/a or for various p (i.e., difxercnt typos of 
machine's) at rated load when a = 1.0. 


Critical dcimplng 

The critical damping conditions for the third 
order system are 

^'2 "" 

™ (B . 9 ) 

pplying these conditions to the system (B,2), two 
relationships 
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3(K-1) (i.Z2)= (B.ii) 

(~^ PK + X-, ) 

'j.rt. olDtrxincd , which on solving, will give the values of 
governor parameters 6 and st the point of critical 
dai.iping. Lqn«(B.10) alono defines the curves corresponding 
to the critical damping point. 
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PART B 

MODAI COKTROI OP 
A POWER S Y B T E frJ 



CHAPTiiR 1 


INTRCLUCTIOE 


In the recent years the optimal linear regulator 
theory for linear, time- invariant systems vith quadratic 
performance indices has been applied to design controllers 
for povor systems for obtaining improved dynamic behaviour 
of these systems [l, 10,11]. Yu et al [l] considered a pov^er 
s^'stem consisting of a synchronous machine connected to an 
infinite bus through a transmission line, along with the 
07 citor-voltage regulator and govern or- hydraulic systems. 

TiiC oi>:',bth order nonlinear mathematical model describing 
this system vras linearized around an operating point to 
obtain a lineer, time- invariant state model of the power 
system, valid for small disturbances. This is represen- 
ted by the standard notation given in equation (l) 

x(t) = A x(t) + B u(t) (l) 

Defining the optimal state regulator problem for this 

system with the quadratic performance index 

‘ 

J = / (S.’ <2 X + u'R u)dt (2) 

o 

where Q and R are symmetric, positive definite matrices. 

Yu et al derived a feedback control law in the form 

u(t) = -R“^ B’ K x(t) (3) 
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Thfci nzn Viain’ matrix Z is the positive definite solution 
to the matrix equrition 

i:' Z + Z A - Z B B' Z + Q s= [Oj (4) 

Thu rociultin^’ closed-loop defined by the equation 

x(t) = (1 - B R-^ B' Z) x(t) (5) 

is a stable system with bettur dynamic response than the 
uncor trolled system e,iv-en by eqn;(l). 

It is to be noted that thure is no apriori knowledge 
of ’h’lut the matrices Q and R should be in order to obtain 
the d'-iiir'eri closed-loop system behaviour. The closed-loop 
system beh^'Viour is judj^od in terms of its transient 
ren]jonse charroct eristics as in classical servo-system 
dusii'ri mc.thods, and this behaviour derends mainly on the 
location of the eie'.envaluus of the closed-loop system. It 
ha.'; thus been necessary for Yu et al to systematically 
change the and R matrices till a satisfactory closed-loop 
behaviour is obtained. It is better to employ design 
tuchnioues which directly yield the required closed-loop 
ei;<cnv.alue locations of the system. This vrould avoid the 
problem of having iteratively to solve a sequence of 
optimal, linear regulator problems as done by Yu et al. 

In contrast to the above design technique, in the 
modal control approach, the desired locations of the 
closed-loop system eigenvalues can be specified and state 
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vr.r.Lra'blc. fc-dback loops ceu be directly designed to shift 
the oigenv.aluus of the open-loop, system to the desired 
locations [2,5 j 4»5]. It is v^ell knewn that it is possible 
to ri,cl;ieve any arbitrarily assigned closed-loop eigenvalue 
locv *;ions (with the constraint that complex eigenvalues 
should occur in conjugite pairs) if the open-loop system 
is controllable [9]- However, in many practical situations, 
one is interested in shifting the positions of only the 
dominant eigenvalues. In these situations it is not 
necor.fc'cry that the sj^stem be completely coiitrollable to 
' chievc this objective; however, it is necessary that 
the odj.,oj:vilueB in que-stion be controllable [5]* lu many 
citup'tions it is possible to aff • ct a shift of the dominant 
ei^,,onvf‘lues with the modal m^..thods to obtain better system 
behaviour with incomploto state feedback [6,7], whereas 
incomplete state feedback cannot be considered directly 
in the earlier technique. Furthermore, tho modal control 
method of design is simple and straightforward . It is 
thus apparent that the modal control approach offers 
f'.dventr'gos over tho method of design, using optimal control 
theory, adopted by Yu et al. 



GHlPTJiR 2 


MODAl- GOl^JTROL APPROACH 

2 .1 iiodal Control I'heor.v and Methods 

Assuming that the eigenvalues of the system 
r‘ TjreOuntod by oqn. (1) are distinct, the time response 
of tiiG systor with zero control inputs in terms of the 
sy: tuhi modes is given by 



■vlicre X and (i - l,2,,,.,n) arc the eigenvalues and 

cif cnvoctors rtispectivoly of the uncontrolled system 
ii'';trix A, c^'s are appropriate constants, and the term 
cxp(a^ t) represents the mode of the system associated 
’■lith eigenvalue X observed that the speed of 

response of the system is determined by the real parts 
of the eigenvalues of the system and the oscillatory 
nature is choracterised by the imaginary parts of the 
ei ^'I’onvalues • So, b5'' shifting both the real and imaginary 
parts of the eigenvalues of the system to appropriate new 
locations in the left half of the complex plane, an 
initially unstable mode can be ma,de stable or the time 
constant of any stable mode can be reduced to the extent 
dv-sirod and improve the overall transient response of the 
system. 
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iiow, for the systoin defined "by eqn.(l )5 it is 
uct-rdbl to shift any specified eigenvalues to any 
-aGr.dyned now locations by linear state feedback (with 
the constraint that complez eigenvalues should occur in 
co'njug.-ito pairs) if the corresponding modes are control- 
I'iblc [ 4 , 5 ]. When a particular eigenvalue, say is 

to b(.' shifted to a. proassigned location, the corresponding 
0 i {'cnvcctor of matriz A’ is to be used as measurement 
otor and the r* suiting output signal is to bo amplified 
by tln-‘ appropriate amount before being fedback. Thus when 
" {:ln) of the oigenvnlues of the system (l) are modified, 

the resulting closed-loop system will have £ new eigen- 
values and £ new corresponding eigenvectors and the rest 
(n-£) eigenvalues and their corresponding eigenvectors 
rt.^'-'in xaichanged [2,3>4»5]» A proof of this is provided 
in AilA.ndix A. 

Thus by using modal control method of approach 
any number of eigenvalues of the system (l) can be 
shifted to some pro-assigned locations without affecting 
th' other eigenvalues in order to stabilize an initially 
unstable system or to improve the transient behaviour of 
rn initially stable system, or both, provided that these 
modes are controllable and all the state variables are 


available for feedback. 
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/;\rious KOdn,! control mi^thods arc diacuss od by 
I'ortor ct al [2,3r4»5] for shifting tho real as well as 
complex oigcnv.vluos of a system repreeented by eqn.(l). 
making use of those methods linear complete state feed- 
i.'-'Ck loops can be designed so that the closed-loop system 
c.-'ii have the preassigned choice of eigenvalues. The 
< Tailed discussion of one of modal control methods .that 
n'’S been used in this work is presented in Appendix B. 

2' • 2 D'..ctor Criterion 

It is essential that the- eitsenvolues of th( closed- 
locp Gyr.tom ; re to be preassigned before the eigenvalues 
of the uncontrolled system are shifted to these new 
3.oc'tiono by using modal control methods. For the pre- 
a.scignracnt of tho new locations of eigenvalues, a sector 
criterion is used which is based on tho fundamental 
concept of damping ratio, ^ and undamped natural frequency, 
w^, ns explained in Figure 1. 

If the dominant eigenvalues of the system that are 

to be shifted happen to be a complex conjugate- pair, th^n 

a. d -.-Sired degree of damping ratio, F, and undamped natural 

Crcqucncy, w^ can bo prescribed and the new locations for 

the conjugate pa.ir can be determined. If another pair of 

corapl^-x conjugate eigenvalues is to be shifted, then ^ and 

w for this pair arc chosen such that the sector defined 
n 

by these K and w^ values is contained within the sector 
defined for the first dominant conjugate pair and the new 
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locf'tiona for the second pair arc placed to the left of 
tkoce preassigned for the first pair. In ca.se a real 
eigenvalue’ is to he shifted, then the new location for 
that r^al eigenvalue can he chosen arbitrarily such 
that it li'. s to the left of the first complex conjugate 
j''‘ir . 

Making use of the above sector criterion procedure 
rif.v locations for those eigenvalues of the uncontrol- 
J wd c^f’stem that are to he shifted, can he preassigned, 

;iO’ 'V'.r, after shifting the uigenvrlue-s to the proscribed 
'.ex" l.oc* tions, if the r'' suiting feodhrick gains arc high, 
tiu-i'i the ciioice of prcassignc'd locations is to he 
u u i t ' ' bly m od i f i od . 

2 . 5 Choice of Input Yariahlos 

When a numb or of eigenvalues of the syrtem (l) 
are to ho siiifted, and if the control variables that 
cfu bo manipulated are many, then the choice of the 
c-'.r.trol input to be used for shifting a particular 
• igonvaluc can be deternined by obsi.rving the terras 

(3 * l,2,..ii, r), r being the number of control 
Inputs and being the eigenvector of matrix A' 
coiTOpponding to the eigenvalue that is to be 

Gi.iftod [6,7]. The inability of a particular input say 
: -tb input, to excite the i-thmfvdc is indicated by the 
I'olationship 



0 


( 7 ) 
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I'iius, it ic! clerxr that there may he more than one control 

input which can excite a particular mode. In such a case, 

a ,;';ood choice can he to select that control input for 
i ^ 

'■•,jich the V h^ term has absolute maximum value among 
all such terms which refer to different control inputs 
tiiat can excite the same mode [6,7]. A different choice 
may lead to a higher feedhack gain requirement. It is 
'jIwO possible; that a particuilar eigenvalue cannot he 
nhil'teid using any of the r available control inputs ; 

;in such a crise the corresponding mode is said to he un- 
cun trollnblc. 

2.4 PoGCihility of Incomplete State Feedhsek 

hpto now it has been assumed that all the state 
vari.'iblos should he available for linear feedhack in 
order to shift the eigenvalues of the uncontrolled 
systCi/ to some specified ne;w locations. But in many 
practical situations all the states cannot he measured, 
and oven if they could he, it maj'" he uneconomical to measure, 
all of them. In such situations it may he possible to use ^ 
incomplete state linear fcedha.ck for shifting the eigen- 
values by modal control methods. A design technique hcis 
been discussed by Davison et al [6,7] for siiifting the 
dominant real eigenvalues of a linecor time- invariant system 
with incomplete state feedback. A similar technique, the 
details of which are given in Appendix C, has been used in 
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tliio work for shifting hoth real and complex eigenvalues 
of the power system under consideration. By using this 
; thod certain unmeasurable state variables as well as 
certain measurable states may be eliminated from the linear 
fetdbr ck and at the same time achieve the proper locationing 
oi tht' oigt.nvalues of the closed-loop system for better 
dynamic behaviour of the system. 

However, it is to be noted tho,t when some of the 
i.i/’ v;nv'"lues ('0 of the open-loop system are shifted to 
pr*. ■•;‘aaignocl locations using incomplete state feedback, 
tin rt suiting closed-loop system may not have the remaining 
(n-’’-) eigonvaluos unchanged unlike in the cose of complete 
ct''ti fcudbnck. The manner in which those remaining 
t igenv- lues get affected hos not yet properly been 
cc t^' blichC'd . Hence, it is necessary to check the cigen- 
V? lue loc''tionB of the closed-loop system after the design 
nd est';blish that the design is satif f-. ctory only when 
tiic .-igcnv'-'luGS which arc not considered for shifting, are 
not i.'Odificd to n large extent. 

2.5 Advnnlages of Modal Control Approach 

Modal control approach can be found to be 
'■■'dv?ntagoous over the method that uses the optimal 
control theory (as used by Yu et al). Computationally 
the modal control method involves mainly determination 
of eigenvalues and eigenvectors, whereas in the latter 
method an .algebraic matrix Riccati equation is to be 
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ooivcc’ using iterative techniques. Purthermore, numher of 
C'OBlgns in an iterative manner are to be done in the latter 
C"..-.e for obtaining a desired degree of transient response 
of the system, whereas in the former case such iterative 
procedures are not necessary. The modal control method is 
simple end otraightforvrard. It gives more- insight into 
t.i.e controlling of the system; that is the control confi- 
gur' tion designed indicates which of the input va.riables 
tn t‘i! nipulate for control and which are the stote variables 
tii.^.t can be used for linear incomplete state feedback. 

'.l'K'.sc advantages become significant when the order of the 
fi;.ot-..' is largo '•■•nd the input variables that can be 
"•■i.i. ulated arc many. Since an implicit control can be 
j/c .'-cisec? over the location of the eigenvalues of the 
closof' -loop system in the modal control method ; any 

dogroc of improvement in the dynamic behaviour 
of the system can be achieved , of course subject to the 
.ngin'.urlng constraints, like gain limitations. 
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CHAPTER 5 


APPPICATIOH OP MOPAL CONTROL APPROACH TQ - THE POl'vER SYSTEM 

3.1 System Model 

The power system whose dynamic hehaviour is to he 
improved by a feedback control law- consists of a synchronous 
genera.tor driven by a hydraulic turbine, having an exciter- 
voltago regulator and a turbine speed-governor as, given- in 
Figures 2,3 and 4. The system was considered- by Yu et alCl] 
for which they obtained a feedback control la,w using the 
linear optimal regulator theory for quadratic performance 
index as defined by eqn. (2). The optimal controller is 
designed with the choice of the following Q e.nd R matrices. 

Q = diag [10 10 1 1 1 1 1 Ij 

r -, 

S = aiag [11] 

Appendix D gives the nonlinea,r equations which 
d-;.scribe the system dynamics, as well as the linearized 
equations derived by linearizing the nonlinea.r equa,tions 
at an operating point. The linearised equations obtahned 
in this manner are 
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Ac- ' rdin^. to tho procoOure outlined in Section 2.3* 
fch* cliOlco of tho control inputs tc shift the five dominant 
. i i/cr.vr-] uos is determined after examining the terms 

1 I 

V hj, (i-i,2,*..,n and j * 1,2). These terms are given 

u 

in TcMo 2. It is found that both the control variables 
u (xcitc all the five modes, but a good choice is that 

li.t ’1 . conirol input is to be used for sliifting the first 
tiiroo dominant oi-onv, clues and the u^ control input is 
tu b' us'-il for shifting the other two eigenvalues. 

If, i i'f- r* lit K-odcl control cinf iguratiens are designed after 
]!rt non ig.ning the now locations for these five eigenvalues. 
I’hM'j ar<; givf.n below. 

control configuration 1 

i-’or the first dominant pair of complex eigenvalues, 
i th the specifications of C =? 0.4 and w^ = 1.0, the new 
lfr''ti<mB e.re determined to be -0.4 +j0.915,in accordance 
■.ith the n-.ctor criterion (Ooction ^2.2 ) . Then for the 
uec'.-nd dominant pair the new locations are taken to be 
-C.96+j0.72r this corresponds to C= 0,8 and w^ = 1,2. 

The ii'jv Icc^^tion for the real eigenvalue is arbitrarily 

nl.QBvn to be -0*6. 

Using the modal control method outlined in 
A'-pf-ndix B, the first three eigenvalues consisting of 
first dominant conjugate pair and one real eigenvalue, 

'■re shifted simultaneously using control. Then the 
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jiKcurui Uominnnt pair of complex eigenvalues is shifted 
u^ control. Ihe resulting feedback gain vectors 

fsrt in Table 4 and the final eigenvalue locations 

riTU t^.i-vvn in Table 5. 

Kodftl control configuration 2 

For this configuration, the new locations for 
th* i Iv*, dtJjT.inBnt eigenvalues of the uncontrolled system 
fu" cf.otj*. n to lie farther to the left in the complex 
I'lrtfn,;, ti.on the corresponding locations chosen for 

1 i ticn 1 in order to see how the feedback gains 
i.’it affect'-d. The preasaigned choice for the new locations 
nlong with the ^ and w^ specifications, is given in 
Table “J* The shifting of the five eigenvalues is then 
carried out in the same order as it is done in the design 
of the confit.uration 1. The resulting feedback gain 
vectors arc given in Table 4 and the closed-loop system 
otgenvaluee in Table 5. 

Modal control configuration 3 

This configuration is designed to show the.t the 
shifting of the first three dominant eigenvalues, very much 
farther to the left in the complex plane, results in high 

feedback gain requirements for the -u control*. The pre- 

s 

assigned choice of new locations is given in Table 3 and 
th© resulting gain vector after shifting the three eigen- 
values is given in Table 4. 
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Tablf 2; 




j=l,2) terms 
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T b, 

mm 
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9 . 5 4 78x1 0~^ 


1 . 0 :; 30,0 
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the Five Dominant Eigenvalues. 


S, Modal Control Configuration 

0 . 1 2 3 4 
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2 

3 -0.60 
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?!o ditl control confifcjiiration 4 

This conf igura'tion. is designed to show that the 
: iii.!. ting of the second dominant pair of complex eigenvalues » 
V(,ry Much farther to the left in the complex plane, 
requii eo high foedbock gain. The preassigned locations are 
( iven in Toble 5 and the r-. suiting gain vectors after 
.^lifting are given in Table 4. 

Mo dal control configuration 5 

This configuration is designed to show that a 
v'Twric choice of control input for shifting a particular 
i.,t.riV"luo may lead to high feedback gain req[uirements . 

* ‘ r shifting the first three dominant eigenvalues using 
of control variable is not a good choice. So, using 
Uy control the three dominant eigenvalues are shifted 
simultaneously. The resulting feedback gain vector is 
|.,ivcn in Table 4. The preassigned locations are the 
snme as specified in configuration 1. 

The transient response characteristics of the 
power system with modal control configurations 1 and 2 
that use complete state feedback are obtained as shown ' 
in Figures 5 to 12. Standard 4-th order Runge-Kutta , 

method is used for solving the differential equations ; 

i 

of the system to obtain the transient responses. The 
initial conditions used are those considered by ! 

Yu ot al [1] which were found from the release of a 
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lond pt.rturbati on The transient response for the 

linoontrollcd system as well as for the controlled system 
\ i 'rh optimal controller (designed by Yu et al.) are also 
obtained for the given initial conditions^ 

3 • 3 I)«.i.^ign, of Modal Control Configurations with 
Incomplete State Feedback . 

It has been decided here also to shift the five 

oo.virifuit cigonv.-jlucs of the uncontrolled system by 

d(... modal control loops with incomplete state 

fo' dback. The preassigned new locations for the five 

doivlnwrit gen values arc considered to be the same as 

tiioru used in designing the modal control configuration 1 

with complete state feedback. The order of shifting 

c 'Hsid' rod is also the same as followed ea.rlier i.e., 

th. first three domin.ant eigenvalues are shifted in the 

bt (.inning using u control and then the other two eigen— 

S 

Viiluea are shifted using u.^ control input. 

According to the procedure outlined in Appendix C, 
first thu average steady state values of the uncontrolled 
system are determined with both the control inputs applied 
in stop function manner [6]. Here equation (C.4) is used 
These average steady state values are given in Table 6. 
Since the first three dominant eigenvalues are to be 
considered first for shifting, the eigenvectors of A' 
corresponding to these three eigenvalues are modified 
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UGin« the £;qn.(C.5). Lxamining the three new eigen- 
vecioro (given in Table 1), it is foimd that the 2nd, 
f.'th, 7th, 4th and 5th elements in each new eigenvector, 

I the smrilli.st, in that order. Since n is already 
ii£j‘.d ar ict.dbock signal in the original system (Figure 4), 
til. otiicr icur states i.e., h, g^, v^ and v^ may be 
(excluded from feedback. 

^iodal control configurations 6,7,8, 9 and 10 with 
iricorijtlute state feedback excluding one (h), two (h and 
gji), three (h, g^ and v^) and(h, g^ and v^) and four 
(h» ££f and v^) states respectively are successfully 
designed which effectively control the first ;hree 
cIc-L.ln-nt modes of the uncontrolled system. The details 
c: the ulgenvaluo location for the resulting closed- 
Icop systems and the gain vectors are given in Tg.bles 8 
to 12. 

Since it has been decided earlier to shift all 
tlu. five- dominant eigenvalues, the second dominant 
c<5njugvte pair is then shifted using u^ control following 
the same procedure as outlined in Appendix C. However, 
ill this case the appropriate closed-loop system resulting 
after shifting the first three eigenvalues with incomplete 
state feedback is to be considered as a starting point. 

Two modal control configurations 6,7 are successfully 
designed with incomplete state feedback dropping one (h) 
and two (hrg^) states respectively, which effectively 
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cohtrol tht fivo ffodes of the uncontrolled system. Design 
of modal control loops excluding three and four states 
irotr. itodback are not successful. The details of the 
i;,ain vectors and the eigenvalue locations for the closed- 
loop ayut^'ms for these controllers 6 to 10 are given 
in Tall 08 8 to 12. 

Tho transient response characteristics for the 
modal controlled system with incomplete state feedback 
dropping two states are obtained, as given in Figures 5 
to 12. 


3.4 Discussion of Eesults 

iiXamino tion of the gain vectors for various modal 
control configurations given in Table 4, indicates that 
in g antral farther the locations to the left in the 
complex plane to which the eigenvalues are to be shifted, 
the more is the feedback gain required. Furthermore, we 
cannot shift the eigenvalues to any large extent be- cause 
of the feedback gain limitations, which is necessary to 
consider from the practical implementation point of view. 
It is also clear from the design of modal control' «onfi- 
turabion that only a proper choice of control input 
for shifting various eigenvalues can give rise to lower 
feedback gain requirements. 

The transient response curves for the modal 
controlled power system with complete state feedback 
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(xij V 'j 11 as incomplGte state feedback, as given in 
Figures 5 to 12, compare well with that of the system 
itii optimal controller designed earlier [l]« In fact 
tiu. response of the modal controlled system with 
coinpl..ti; state feudbo.ck is better than that of the 
optimal controlled system designed by Yu et al [l]. 
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S Q 

iftblc 6i X — Average Steady State of the System 

with "both Inputs given in Step 
Function Manner. 


S.iio. 

State 

Value 

1 

6 

1.1605 

2 

n 

0.0 

3 

’^f 

-18.4720 

4 

^f 

-2.7199 


"^s 

1.8248 


g 

-0.. 8117 

7 

^f 

0.0365 

8 

h 

0.0010 


Tabic 7 j Modified Eigenvectors for the First Three 
Dominant Eigenvalues. 


S.No. 


v^ 

1 

-0.5981+j0.1208 

-0.0292 

2 

0*0 ® 

0.0 

3 

1.1730+j0.2082 

-0.0188 

4 

-0.0006+;j0.0297 

-0.0048 

5 

-0. 0311+ jO. 0095 

-0.0202 

6 

-0.5548+10.1985 

0.8075 

7 

0.0048ij 0.0006 

0.0055 

8 

0. 0010+ j 0.0 

-0.00002 
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Tabic 9.* Xodal Controller with Inconpl.tc Jt-t-i 
excluding ’h*, 'g„' states (ccriif.urati 
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12.8758 
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e- With modal controller config.l 
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( excluding h,gx states) 













CHAPTHiR 4 


OORCLUSIOR 


Modal control method of approach is used in this 
work for designing feedback controllers for a power system 
in order to improve the dynamic behaviour of the system. 
Modal controllers are designed with complete state feedback 
as well as incomplete state feedback, which effectively 
ir.r)rovc the rcsy'onso of the system. The advantages of the 
rcdal control method over the earlier method which uses 
optimal control theory arc discussed. 

borne of the possible extensions of the present 
r«, search are listed below: 

(i) The power system model considered in this work 
is of 8-th order consisting of third order machine, second 
order exciter-voltage regulator and fourth order prime-movei 
with speed governor. More detailed dynamics of the machine 
as well as prime— mover can be included in the model for 
modal control study. 

(ii) In this work the modal control approach is 
.-nployed to control a single machine power system. The 
same technique may be applied to control a multimachine 
power system. 
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APPENDIX A 


Por a linear time-invariant system represented "by 
the equation (1) 

x(t) = A x(t) 4- B u(t) , 

let US suppose that the eigenvalues of the open-loop 
system arc distinct. Let D^, (i=l, 2, . . . ,n) are the 

eigenvectors of matrices A and A' respectively, and 
difigonal mo,trix A is formed by the n eigenvalues 
(^ 1 > ^ 2 '***^ system,. Then 


A = U A V 

V A U » A 

A ^ i 2^ 

A’Y^ ^ (i = l,2,..-.,n) 


(A.l) 


(A. 2) 


The two sots of oigxnvectors are orthogonal, i.e., 

= 0 (i ^ j; i,-j=l,2, . . . ,n) (A. 3a) 

ar.c' can bo normalised so that 

1 (i = l,2,...,n) (A. 3^5) 

Suppose the first r eigenvalues of the system (l) 
arc shifted using r control variables. Then the first 
r eigenvectors of A’ are to be used as measurement 



81 


vectors and the control law is given by 

^i ” ^i X (i = 1,2, . . , ,r) (A. 4) 

where are the controller gains. The closed-loop 

system is given by 

/V A 

i(t) = (a + B k V’ ) x(t) 

= A^ x(t) (A. 5) 

A 

whore k is a (rxr) diagonal matrix formed by the 

A 

C03'j troll or gains and V is the (nxr) matrix given by 
V = [ ] 

In view of eqns. (A. 2) and (A. 3) 

= A 2.^ + B k V* 

*= k^ (i = 1,2, . . . ,r) 

= A 

= (i = r+1, r4-2, .. . ,n) (A. 6) 

It can be deduced from the eqn.(A.6) that the last (n-r) 
eiponvalues and eigenvectors of the closed -loop system 
iiTitrix At|_ arc the same as those for the open-loop system 
I'T'.trix A, and only the first r eigenvalues and eigenvectors 
get modified. 
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appendix b 

MODAX CONTROL JV-LETHOD [4,5] 


Por the linear time— inve riant system represented 
by eqn, (l) 

x(t) = A x(t) + B u(t) , 

lot -us suppose that i of the n eigenvalues either real or 
complex, or both arc to be shifted simiiLtaneously with the 
rt.Gtriction that complex eigenvalues occur in conjugate 
pairs. This can be achieved by using a single control 
input [4,5]. Let u^ control input is selected from 
among the r available control variables (in accordance with 
the procedure given in Section 2.3). Then the closed-loop 
system is 

x(t) = A x(t) + b^ u^ (E.l) 


Let i feedback signals arc generated from the 
.’.le.'isurement of all the states using eigenvectors 

jL 

(i =5 l,2,...,il) of A' as measurement vectors and are 
-mplified by a respective controller gains k^, (i=l,2 , . . . , 
Then the control law is given by the sum of these £ signals 



(E.2) 


The closed-loop system is 
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a 

x(t) = (a + 1^2 ) ^(■t) 

^ A*! ) 


(B.3) 


In view of the orthoganality conditions of the 
eigenvectors as given in eqns.(A»3)» it follows that 


a 

= A U^" + h^ kj, 

= ( j=a+l,a+2, . .. ,n) 


= Xi U" + h^ k. 


(i= 1,2, . . . ,£) (B.4) 


Proni uqn.(B.4) it can he deduced that when a eigenvalues 
of the uncontrolled system are shifted, the resnlting 
closed-loop systci' will have a modified eigenvalues and 
xheir corresponding eigenvectors and the remaining (n-a ) 
ejgenva.lues £ind eigenvectors will he \inchanged from those 
of the open-loop syster.. 

In order to calculate the controller gains 
kj^ (i = l,2,...,a), let the new eigenvalues and eigen- 
vectors he Pi and ( i = 1,2,..., a). Then 

(B.5) 

(i=l,2, .. . ,a) (B.6) 


= X^ 


lot 


rind 


n 


# = I li-i 2^ 
2=1 


n 


= I P-5 2 


2 


(B.7) 
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’’horw and Pj^s are some constants. In view of the 

eqns.(ii.5) and (B.7)» the p. constants are given hy 


(3=1,2,..., Jl) 


(B.8) 


Substituting from cqns.(B.3), (B.6) into eqn.(B.5) yields 

('•- + >^1 1, i"’ ) 2 ’ = Pi hxi 2 ^ 


( i— 1 ) 2.f r • * f z) 


(B.9) 


Since 




I P-li = f 4.1 \ 2^ (i=1.2 *) (B-10) 

3=1 j=l J 

hleo, in view of eqn.(A.5) and (B.8) 

fli 1 li 

3=1 


|/.i 2 ^ L 


><1 U 

=1 0 f=l 

( i=l , 2 , . . . , if-) (B . 11 ) 

Making use of eqns.(B.lO) and (B.ll), eqn.(B.9) can he 

) f 

written as 


n . n ^ 

I + I Pj 2.^ I = 0 I Q.- - 

3=1 ^ 3=1 3 f=l ^ ^ 1 jii 13 - 


(i = 1,2,...,^) (3.12) 
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Eqn. (B.12) is a vector equation in and is equivalent 

to the (nx®') scalar equations 

I 

<Pi - iij - Pj = ° 


i=l 0==:l,2,t..,n 


(B.13) 


For a iiiven value cf i, the first i of these (B.I 3 ) 
e 'a' t.i ons , i.e.^ those corresponding to the eigenvalues 
to be shifted, can be written in the matrix form 

% = 0 (B.14a) 

wLi rc 


H. 




(B.14b) 




~ ( i > 0 » f~l » 2 ) 


(B.14c) 


(B.14d) 


and 6.,. the Kronecktr delta. Since q. ^ 0 it 

j 1 ^ 

follows from {B.14a) that 


j 


0 


(B.I5) 


In view of equations (B.14b) and (B.14c), this has the 
explicit form 

; '■^1^2 *'* ] 


"i 2' -^2^2 


‘V% 




(B.li) 
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laxpansion of the dotorminant in oq,n, (B.16) yields 
Z a ^ 

® ^ ^f ^f ^ (i=®li 2 , . f . ) (B.17) 

d=l f=l j=l 

which implies 

^ P ^ 

j _J — J — = 1 (i=l,2,.*.,£) (B.18) 

0=1 (p j ) 

provicod that 

y 0 (i,o = 1,2,..., 

Bqn. (B .18) can he solved for to give 


k 


d 


z 

TT (p^» 
f=l ^ 



Z 

IT 

f=l 

f^j 


(x,-V 


( 0 = 1 » 2 ,. 


il) 


(B.19) 


vrhich indicates that the controller gains k^- can he 
rt.j;lisc-d only when 

^ (i,0 = 1,2, ...,!l ) (B.20) 

r..nc' 

0 (i = 1, 2, . . . , ^- ) (B.2l) 


The inequality in eqn.(B.20) implies that the eigenvalues 
to he shifted should he distinct and the inequality given 
hy eqn.(B.2l) implies that the first it modes of the un- 
controlled system shotild he oontrollahle. 
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Th'. control law obta,ined by substituting eqn, (B,19) 
1 -> nqn,(B.2) has the form 


u. 


I 

I I 

i=l 




Pj^ ^IT ^ ^ 

f^i 


ft 

I X ] 


= G’ X (B.22) 

t'(. £ iH thu gain vector the elements of which are 
. i ^ *.r. by 

2 , 

G = I I t' (B« 23 ) 

i=i p n (x^-x ) 

^ f=i ^ ^ 

f^i 


It is to be noted that in case of real systems 
thn control law d>;fined by eqn. (B.22) will be real even 
in the erso of complex eigenvalues. This follows from 
th(, fact that if X^ ~ x^* , then = V^* and Pj_ = p^* , 
together with the fact that the required new eigenvalues 
^ ^ rill cither occur in conjugate pairs or be real, 
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APPElvDIX C 

DfeSIGB TECH I jlQUJ:; FOR niCOMLETE STATE 
FBEDBAGL OOKTROLLBR 


let us s^ippose that I eigenvalues of a linear time- 
invoriant system as represented hy eqn.(l) 
x(t) = A x(t) + B u(t) 


ru’' to he sidfted to preassigned locations with incomplete 
otto feedback using a single control variable, the choice 
M’ •. :ach is i.'.ade according to the procedure outlined in 
c-tion 2.5. Then Si measurements of the state variables 
■ vfi t' bo mode using the Si corresponding eigenvectors of 
A* ns mofisuroment v -ctors; and these i signals are given 

by 

T,^(t) = v^’ x(t) 

(i = ) 


The quantity gives a measure of how much of the i-th 
mode is present in the system at any instant of time. 
This follows since if x(t) is written in the form 


? i 

x{t) = I E 

~ i=l 


uhen 


r)^(t) = t £(■*) = 


(C.2) 

(C.5) 
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! t. n.fiy h'lppon that under normal operating conditions 

of nrc much smaller than the others, then 
'.If.. cun brj omitted from measurement without affecting 
'/."Jui. cf significantly, V/ith this basis of 
'■'! '• -liii utic'n of the modes, the choice of best state vari- 
’ i*' ir incomplete state feedback can be determined[6,7] • 

T tr;C(,tlurc of design is as follows: 

..t: X ** »^otormino the a eigenvalues and the corresponding 

‘ i, • hv* ctors of the transpose of the system matrix. 

U'-t rmino' the I®® - the average steady state 
V" I for the system state vector when the input forcing 

roncti'-rta are applied in a step function manner [6]. These 


vrj uer will normally correspond to typical operating 


’/••'u'.o of states, at least when there is no control applied. 



r 

I 

i=l 




* h 
i — i 


(C.4) 


o tMi 3 : ilcdify the i eigenvectors using the steady state 
values. The elements of these new eigenvectors are 


v^' » x!®v^ (J-l,2,...,n; i«l,2, . . . , (C.5) 

J j 3 

Stop 4 : ii)xamine each of the new eigenvectors for elements 
whoso absolute value is smaller than that of others by an 
order of 5 to 10. Select those states for which the 
corresponding elcm.^nts of the new eigenvectors arc the 
smallest i and at the same time these states should be 
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C'.'.rfnron t * r.ll the now eigenvectors. For example, suppose 
p,q 'ind r tltments arc the smallest in the first new eifeen~ 
V'-Ct r V find in the subsequent eigenvectors only p and q 
• r. t!. nraf'lli.st elements, then only p and q states are to 

1/ i M , c t ^ 4 ♦ 

Thonc st^'to v-'riablcs selected in the- above, manner 
C'l; b« .mi tied from feedback and thus we may achieve 

.otnto feedback control of the system. At this 
t,l-n V(. hope for eliminating unmeasurable states as 

v ».o £3' me vf the mea-surablo states. 
t.<t-.s 5 i Obtain the control law assuming complete state 
fe' (ibfick ucing the modal method outlined in Appendix B. 
utep 6 ; heglcct the states as determined in Step 4 from 
fi.edbnck and obtain the controller with incomplete state 
Ic dback. 

htop 7 t Obtain the new steady state values applying this 
c-ntrol lav ajid chock how close those true steady state 
values Jirt tc. the values calculated using the eqn. (G.4)* 

If there is a substantial differ once in the relative 
ratir- of various elements of then tha steps 5 to 

6 art. t> bo ropfcuted* 

Stop 6 : Finally obtain the eigenvalues of the resulting 
closed-locp system with incomplete state feedback and check 
how the eigenvalues which are not considered for shifting, 
aro affected. case they are affected to a great extent, 
then the dosisu may not give satisfactory results. 
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Gti,.,p 9 » T/^hcn a niiinbor of eigenvalues of the system arc 
to Lo shifted using different control inputs, this can he 
d'.no sequentially repeating the above 1 to 8 steps, each 
tinu rjhifting the appropriate eigenvalues with the corres- 
loiiding choice of control input. Hox^ever, it should be 
noted that in this sequential procedure same state 
vir:r i'lcs are to be excluded from feedback at each sta.ge. 
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API^El^DIZ B 


Nonlinear equations of the power system 

The power system considered by Yu et al (1970) 
is represented by a system of eight nonlinear equations 
which are given below together with the linearized 
equations and the numerical values assumed. 


Synchro n ous machine equations 


w 


w 


- D w - p;) 


wh' ro 




% ’-^f 

^d \lo 






“V „ cos 6 


^d 


(D.l) 

(D.2) 

(D.5) 


P! 

1 


=! g + 1.5h 


v^ sin 6 
^do 


(xl - X ) 

^ d q. 

2Zi X 
d q 


Vq sin 26 


- P 


o 


Lxcitcr-voltago regulator equations 


Vh 


V 


V ^ u 

£ — £ (v’' + V ) 

X X ^ t S 

e 0 

V u 

_ _s + Is 

T . 1 


(P.4) 

(D.5) 


s 


s 


Govemor-hydraulic system equations 


* cr g 1 /■ w \ 

e = - T " ■V* “ 


(D.6) 
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+ 


a 

2a g 
6 


_a 

'a 

2 


u. 


•W 

w. 


Ef . 


2h 


w 


(D.7) 

(D.S) 


The nonlinear equations (D.l) to (D.8) of the power system 
are linoerised around an operating point and the real time 
is scaled for the convenience of computation, which is 
given by 


i: = p t 

where p = 7. 308. 


(D.9) 


Liner'rlsed equations of the power syste m 

The lincorisod system equations with scaled 
time are given by 


W 

x(i:) = A x{i) + B u(x) 


(D*10) 


wh(-re 


0 


0 0 0 0 0 0 


-D 


^21 pM ^23 


p^M 


0 


1>5 ■ 
P^M 


ii= 


a^l 0 a^^ p 


0 0 0 0 


0 


w t 
o g 


0 


0 0 


0 


0 


0 


0 0 0 0 


^41 ° ®'43 Pt 


zi 

pT 


0 0 


0 0 


0 0 


0 0 

-a -1 

Bt B% 

^ S ^ E 


0 0 


hi 

2<J ^ 


0 

0 

0 

0 

-2 



94 


°21 


cos 6 ^ ^ (xi - X ) 

[_° O_jo ^ cu_ 2^5 ] 




d do 




a 


23 


v^ sin 6 
o 0 

^d •'do 


a 


31 


(Xd ” xp Vq sin 6 ^ 

B Z’ 

P "cl 


a 


33 


P I’ t:’ 
d do 


a 


41 


\x V X V, cos 6 x' V sin 6 
e 0 j- q do ^ , o d go o 




o"to 


^d 




"43 


X Y 


_go 


P V, ^ XJ "T^i 
*^ 0 X 0 "d do 


B 


0000100 0 


0 0 0 O'O 01 0 < 


x(x) = [6 n 4 'f Vg g g^ li]‘ 


and 


u 


("c) = [u„ 


gj 


where 

A ^‘B 


S- “ P T 


s 


A ""2 

"^g p Y 


Besides the other variables are defined by the following 
equations 
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=: (I-XB) + X 

X = (1-xB) x„ + X 

q. q 

= (1-xB) x^ + X 






+ 



V‘ 

q. 


d 


V' 

q 


% ^0 °°^ ^0 

X ''*„ x' sin 6^ 
Ido 0 


^d"do 


(D.ll) 


J'nrther details of the power system are contained in 
the paper hy Yu et al (1970), 

Nu merical values for different parameters 
Transmission lines x = 0.7417> B = '<'“.1339 
Synchronous machines x^ = 1.0, x^ = 0.27, = 0,6 

t^^= 9.0, M=0. 02122, B = 0.00537 
hxciter-voltage regulators p,^ = 10.0, r^^l.C, '1:^=0. 5 
Governor- hydraulic systems a=0.045, 1:^=0. 1, Tg^=0,01, t^=1.6. 

and values are not needed since they are 
included in u^ and u^ respectively. 

For an initial load and a terminal voltage of the 
synchronous machine, 

Pq = 0,735, Qq = 0.034, = 1.05 

the initial currents, voltages, flux linkages and torque 
angle are 

^do "" 0*286, iqo = 0.64, v^^ = 0.384, v^^ = 0.977, 
v^ = 1.263, v^ = 1.058, = 9.491, 6^ = 0.887. 
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